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Abstract 



> 

\D . Studies of new hyperbolic systems for the Einstein evolution equa- 

tions show that the "slicing density" a(t,x) can be freely specified while 
the lapse N = ag 1 ^ 2 cannot. Implementation of this small change in 
the Arnowitt-Deser-Misner action principle leads to canonical equations 
that agree with the Einstein equations whether or not the constraints 
are satisfied. The constraint functions, independently of their values, 
then propagate according to a first order symmetric hyperbolic system 
whose characteristic cone is the light cone. This result follows from the 
?— I ■ twice-contracted Bianchi identity and constitutes the central content of 

the constraint "algebra" in the canonical formalism. 

In this papeiQ I will show one way to follow a path to some small, but striking 
. improvements in the Arnowitt-Deser-Misner (ADM) canonical form of general 

relativity, jl], || Arlen Anderson, collaborator in the primary reference H on this 
subject, and I followed a different track in the presentation given elsewhere. || 
Several aspects of the ADM canonical form of general relativity jl], |2| have 
been improved recently. || Slight modifications give rise to possibly significant 
changes of perspective. The modifications were suggested by the structure of 
the evolution equations written in new, manifestly hyperbolic, forms, not di- 
rectly relevant to the present work, in which only physical characteristic speeds 
appear. [|| |[ ||, fj], |[ ||, |l(], |ll) One may see some aspects of those results as 
flowing to, or from, the action principle. 

Here I focus on the consequences of a direct demonstration that the con- 
straint functions, whether they vanish or not, always propagate as dictated by 
a first order symmetric hyperbolic system whose only characteristic cone is the 
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full physical light cone. This system is in fact just a new reading of the twice- 
contracted Bianchi identities. This reading suggests the form that the canonical 
momentum evolution equations ought to take. This observation, in turn, causes 
one to change the treatment of the lapse function N; and this change leads to a 
small alteration in the action principle. One's view of general relativity is then 
changed by noting that parameter time differentiation of the canonical variables 
now holds correctly throughout the entire phase space, among other things. 

The definition of the Einstein tensor = R^ v — \Rg^ v in terms of the 
Ricci tensor leads to the observation that 2G[j = Rq — R% and to the 
important identity 

Gij + gijG = Rtj — QijRk ■ (1) 

The vanishing of the right hand side does not depend on any constraints and is 
equivalent to Rij — 0. Evidently, Rij — is also equivalent to Gij = —gijG^. 
Therefore, while i? M „ = and G^ — are equivalent, Rij — and Gij = are 
not equivalent equations of motion unless the Hamiltonian constraint 

H = g x /' 2 C = 2g 1 ' 2 G° = g 1/2 {K %] K l] - H 2 - R) = (2) 

holds exactly. (Here, H — Kj denotes the trace of the extrinsic curvature 

of a spacelike "time slice" with metric g^, density g 1 ^ 2 = (det^) 1 / 2 , scalar 

curvature R, and spatial covariant derivative Vi.) 

These observations lead immediately to a transparent form of the twice- 
contracted Bianchi identities V^G^ = 0. Namely, 

VpG$ = VoGjj + V 3 G Q = , (3) 
V/jGj = V G° - VjG° + Vi [R) - 5)Rl] = . (4) 

Recall that the momentum constraints are given by 

Hi = g 1 ^ = 2g 1 / 2 NR° = 2g 1 / 2 NG° l = . (5) 

Let us write the Bianchi identities (^) , (||) in 3 + 1 form using the spacetime 
cobasis indicated by parentheses in the spacetime metric 

ds 2 = -N 2 {dt) 2 + gij (dx i + (3 l dt) (dx j + (3 J dt) , (6) 

where N is the lapse function and (3 l is the shift vector. The dual vector basis is 
da = dt — ftdj, with dj = d/dx^ and dt = d/dt. We define our time derivative 
operator on time-dependent spatial tensors by <9o = (over-dot) = dt — £ p , with 
£p the spatial Lie derivative along the shift vector. We also use the spatial 

covariant derivative Vj. 

From the cobasis {6° = dt, 9 l = dx l + /3 l dt} and d9 a = -\C a ^ A we 
find C l oj = —C i jo = djfl 1 and all other structure coefficients vanish. This gives 
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the connection coefficients 



7 jk ~ 1 3 k ~ 1 jk 

l\j = -NK\j 

7% = -NKij + di? 

7 l 00 = WW (7) 

7°oi = 7 ,o = ^log7V 
7 ° 00 - a lo g 7V, 

where T denotes a Christoffel symbol. Our derivative convention is V 0,073 = 
d a <Jp - v P l p p a - 

Now we can go through a straightforward process of writing (|3|) and (|4j) in 
terms of Bq ( ) = ( ) and Vj. For the unweighted constraint functions C and Cj, 
the Bianchi identities then become 

C-NV j Cj = 2(CjV j N + NHC - NKij[TZij}) , (8) 

Cj-NVfi = 2fcV j N+^NHC j -V i (N[R ij ])j , (9) 

where TZij = — 9ijRk- Similar results hold for H — g x ^ 2 C and Tij = g x ^ 2 Ci. 
(Related formulas were found by Choquet-Bruhat and Noutcheguemc [O in 
studying the evolution of matter sources p 00 , p 0z , where p@ a = T^ a — ~gP°^Tff .) 

The system (g), @ is first order symmetric ( "symmetrizable" ) hyperbolic 
with characteristic fields consisting of linear combinations of C and Ci propagat- 
ing along the light cone. Hence, if the equations of motion TZij = (or Rij = 0) 
hold, C and Ci remain zero in the physical domain of dependence associated 
with the region of the initial data surface on which they were initially satisfied. 
If none of C, Cj, TZij is exactly zero, as is inevitably the case in approximations 
and numerical applications, then @ show how errors in C and Ci are driven 
along the light cone by terms linear in these errors and by Rij — gijR\- (It is 
straightforward to add a cosmological constant and a matter source T Q/3 , such 
that VpT^ = 0, to this analysis.) 

To draw some lessons for the canonical formalism, let us first express the 
3 + 1 evolution equations in their standard geometrical form (with zero shift, 
Ref. 13; arbitrary lapse and shift, Ref. 14; spacetime perspective, Ref. 15): 

9ij = -2NKij , (10) 
Kij = N (-R^ + Rij + HK t j - K ik Kj - N^VidjN) . (11) 



A brief look at (^Tj) shows that the term TZij = (Rij — gij^k) that appears 
in the hyperbolic form (@), (ph of the Bianchi identities is made-to-order to 
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produce an equation of motion for the ADM canonical momentum 



TT ij = g 1/2 (Hg lJ - K ij ) (12) 
that contains no constraints. Indeed, using ( |l0| ) and ( |Tl] ) , we obtain the identity 

7T iJ = Ng 1 ' 2 (Rg zj - R ij ) - Ng- 1 ' 2 - rnr 1 ^ 

+g 1/2 (V l V J N - g ll V k V k N) + Ng 1 ' 2 [Jl lJ ] . (13) 

From (^J), we have 

9ij = Ng- 1 ' 2 (27r y - - n 9ij ) . (14) 



If we now compute the time derivatives of C (g, ir) and C4 (g, 7r) using (13), (|14|), 
we obtain, of course, the Bianchi identities in the form (||), ([)[). This merely 
stresses that the Bianchi identities, properly construed, are just the hyperbolic 
evolution equations for those phase space functions whose vanishing yields the 
"constraint hypersurface" in phase space, and in particular we see that (|J), (||) 
hold on or off the constraint hypersurface. 

We have now reached a crucial point in the development. If the canonical 
equation for 7T y is dictated by the vanishing of the spatial part of the Einstein 
tensor, G y = 0, as in the ADM analysis, |jj then the identities ([!]) and ( ft3| ) 
show that the constraint term (1/2) g Jl ^g 1 ' 2 C = (1/2) g^Ti remains in the -h l i 
equation, restricting its validity to the subspace of phase space on which the 
constraints are satisfied (i.e., where the constraint functions vanish). Further- 
more, if we substitute G u = back into the Bianchi identities (||) and (Q), then 
the hyperbolicity and well-posedness of the constraint evolution are lost. (The 
latter has been observed by Frittelli using different methods. [^6j She also found 
that = gives well-posed evolution.) 

However, though the ADM derivation of the tt 4 - 7 equation, found by varying 
gij in their canonical action (167rG = c = 1) 

S [g, 7r; N,(3,) = J d A x (^ 9ij - NU) , (15) 

with N(t,x) and f3 l (t,x) as undetermined multipliers, is of course perfectly 
correct, another point of view is possible. (We are ignoring boundary terms, 
which are not of interest here, and we note that the momentum constraint term 
—{PHi (Hi = g x / 2 Ci) is contained in tt^ g^ (( ) = do) upon integration by 
parts.) The other point of view arises in our work on hyperbolic systems with 
only physical characteristics. |Q, |], p[ M, ra, @, |h| [llj There, one is never allowed 
to specify N(t, x) freely! Rather one must use in essence Choquet-Bruhat's 
"algebraic gauge," which asserts that the weight-minus-one lapse (the "slicing 
density") a = Ng^ 1 / 2 = a(t,x) > is freely specifiable while N is not. ]Tt[ | 



(The slicing density a has been used prominently in the action by Teitelboim [118 
and Ashtekar |19l E0| for other purposes 
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Indeed, if one computes do log a = f(t, x) from a given a(t, x), then 

g 1/2 d a = d N + N 2 H = Nf , (16) 

which is the equation of harmonic time slicing with f(t,x) = <9ologa acting 
as a "gauge source." [^l] In this sense the undetermined multiplier a(t, x), like 
l&it, x), cannot affect the issue of hyperbolicity. §, |, | |, @, |, § |g jij The 
"harmonic time slicing gauge" is not a gauge at all; it is simply the equation of 
motion of the lapse function (with a "gauge source" ) once one recognizes a as 
the true undetermined multiplier. 

Combining ( |l6| ) and the equation for H obtained from (|ll|), one obtains a 
quasi- linear wave equation (hyperbolic with characteristic speed c = l!) for N . 
Therefore, every foliation of a globally hyperbolic spacetime by regular time 
slices is given by some initial time slice, and the solution N > of a wave 
equation, for some a(t, x) > 0. No slices are "lost" by changing from N to a, so 
the latter is as "good" as the original N. (Degenerate cases can be handled as 
well with g 1 / 2 and a as with N.) We conclude that N is a dynamical variable 
(as did Ashtekar for other reasons) which determines the proper time 

NSt between slices t = t' and t = t' + St. It is determined from a(t, x) and the 
solution g 1 ! 2 of the constraint equations. |l5|. ^2|, [231 ^i) From this perspective, 
the Hamiltonian constraint plays its familiar role as an initial-value constraint 
that determines g 1 ^ 2 given the other free data. |l5| The important point is that 
ri = does not determine the time but does fix the rate of proper time r with 
respect to parameter time t: dr/dt = g x ^ 2 a along the normal do- 

Motivated by the above findings and those with respect to the ?r y equation 
and the well-posedness of the Bianchi equations as a causal hyperbolic system, 
we alter the undetermined multiplier N in the ADM action principle to a. Thus 
the Hamiltonian function becomes 

H = g l ' 2 H = 7r tf 7T« - \i^ 2 - gR , (17) 

where TL is a rational function of the metric of scalar weight plus two. The 
ADM-plus-Teitelboim-plus- Ashtekar action becomes 

S [g, tt; a,/3) = J d A x ( -••'//„ - afi) . (18) 

The modified action principle for the canonical equations that we propose is to 
vary 7r u and gy, with a(t,x) and ^{tjx) as undetermined multipliers. From 

5H = (27Tij - gijir) Sw ij + {^^^l - jnr« + gR lJ - gg lj R) S gij 

-g (V'V'%,, - g ij V k V k S gij ) (19) 
we obtain the canonical equations 
57~t 

9io = a ^IJ = a ( 27 % - *9ii) = -2NK i: j , (20) 
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sn 

-a- — = -ag [W J - Hg UJ ) - a I Xir"""?^ 



(R lj - Rg lj ) - a (2n ik ir{ 



figij 

+g(V l V j a-g i] V k V k a) . (21) 



Equation fl2J) for tt ij is the identity (|13|) with TV 3 = 0, that is, R» = 0. 
Thus, (|2l]) is a "strong" equation unlike its ADM counterpart which requires in 
addition the strict validity of a constraint: Ti, = 0. 

In the present formulation, the canonical equations of motion hold every- 
where on phase space with any parameter time i, a necessary condition for the 
issue of "constraint evolution" to be discussed properly at all in the Hamiltonian 
framework. Further, there seems to be no need to vary a and (3 l in a spacetime 
action volume integral to enforce the constraints for all time, because the con- 
straints are dynamically determined to hold in the appropriate physical domain 
of dependence if they hold initially, as shown below in (|2^), (|27|). Perhaps a 
and (3 % should be varied on an initial slice only. 

If we define the integrated or "smeared" Hamiltonian constraint as 

H a = [ d 3 x'a(t,x')H , (22) 



the equation of motion for a general functional F[g, 7r; t, x) anywhere on the 
phase space is 

F[ g ,7r;t,x) = -{H a ,F}+d Q F , (23) 

where ( ) denotes our total time derivative and do is a "partial" derivative of 
the form dt — £ ^ acting only on explicit spacetime dependence. The Poisson 
bracket is 

\FG\= f d 3 x ( SF 6G — SF \ (24) 

J X \5 9ij {t,x)5nv(t,x) Sg l3 {t,x)STr^(t,x)J ' { ' 

It is clear that the (g, 7r) equations come fro m (p3f ) applied to the canonical 
variables. The harmonic time slicing equation (|16|) results from application of 
(23) to N, and the wave equation for N comes from a repeated application of 
©to©. 

Time evolution is generated by the Hamiltonian vector field 

Xn a = ( d3x \ a ( 27T v ~ - [ag(R 13 - Rg %3 ) 

J ( ogij 

+ a(2TT lk ir J k - mr ij ) - ,g(V l V J a - 5 y V fe V fe a)]^-| . (25) 



Because is does not contain any explicit constraint dependence, (25) is a valid 
time evolution operator on the entire phase space. 
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The product rule applicable to do and the Poisson bracket shows that the 
evolution equations (^), (^) for the constraints can be written as 

d H = - \n ai n} = agg ij diHj + 2gg tj H i V j a , (26) 
d Hj = - \u a ,Uj } = adjU + ZHd s a , (27) 

where VjCe = dja + ag~ 1 ^ 2 dig 1 ^ 2 . Therefore, the Poisson brackets of the 
smeared with the unsmeared constraints are well-posed evolution equations 
for the constraints. These equations are, of course, just the twice-contracted 
Bianchi identities when IZij — or Rij = 0. 

The consistency of time evolution with respect to different choices of a is 
checked by using the Jacobi identity, 

|w ai , |7i Q2 ,F|| - j7Y Q2 , |7Y q , 1 ,f|| = | , 7^ Q2 1 , f| , (28) 



and 



{H ai ,n a2 } = - J d 3 xgg 11 {a 1 d l a 2 - a 2 d l a 1 )H J . (29) 



The metric dependence of (29) shows that the difference between evolution with 
o?2 followed by ai, and vice versa, is a spatial diffeomorphism when TLi = or 
when 6F/6ir lj = 0. 

The results of the previous paragraphs shed new light on the Dirac "algebra" 
of constraints (cf. Ref 25). It is well known that the Dirac algebra is not the 
spacetime diffeomorphism algebra. This can be seen from the fact that while 
the action (^) is invariant under transformations generated by the constraint 
functions even when they do not vanish, |26| the equations of motion that follow 
from this action when the constraint functions do not vanish are Rij = 0. These 
equations are preserved by spatial diffeomorphisms and time translations along 
their flow in phase space, but a general spacetime diffeomorphism applied to 
Rij — mixes in the constraints. A comparison of (@), (JsJ) with (|26]), ( p7| ) shows 
this effect: the Bianchi identities are spacetime diffeomorphism covariant while 
(p6|), ( p7| ) are not. The equations (||), (||) and (p6|), ( ^7j ) differ precisely by terms 
proportional to IZij and V 1 ??.^-. 

A second important view of the Dirac algebra results from the direct and 
beautiful dynamical meaning of its once-smeared form. Equations (|2^) and ( p7| ) 
express consistency of the constraints as a well posed initial-value problem. If 
the constraint functions vanish initially, they continue to do so by evolution into 
the domain of dependence corresponding to the region of the initial time slice 
on which they initially vanished. This mechanism follows from the dual role of 
Ti. as a constraint and as part of the generator of time translations of functionals 
of the canonical variables anywhere on the phase space. 

Let us take note that the Hamiltonian constraint per se does not express 
the dynamics of the theory; the equation of dynamics is (p3). In its "altered" 
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role, the Hamiltonian constraint function simply vanishes as an initial value 
condition, from which 

ff Va i s 

determined as in the intial value problem. pl5[ 
Then TV can be constructed from a. The Hamiltonian constraint, once solved, 
remains so according to the rigorous result embodied in (|26|), (p7j). 

The application of these ideas to canonical quantum gravity will appear 
elsewhere. [ p7| 
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